// 数值不稳定算法
#include "stdio.h"
#include "math.h"

int main( )
{
	double y_0 = log(6.0 / 5.0),y_1;
	int n = 1;
	printf("y[0] = %-20f",y_0);
	while(1){
		y_1 = 1.0 / n -5 * y_0;
		printf("y[%d] = %-20f ",n,y_1);
		if(n>=40) break;
		y_0 = y_1;
		n++;
		if(n % 2 == 0) printf("\n");
	}
	return 0;
}